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Abstract 

We show that the autocorrelation of quantum spectra of an 
open chaotic system is well described by the classical Ruelle- 
Pollicott resonances of the associated chaotic strange repeller. 
This correspondence is demonstrated utilizing microwave ex- 
periments on 2-D n-disk billiard geometries, by determina- 
tion of the wave-vector autocorrelation C(k) from the exper- 
imental quantum spectra S2i(k). The correspondence is also 
established via "numerical experiments" that simulate S%i(k) 
and C(k) using periodic orbit calculations of the quantum 
and classical resonances. Semiclassical arguments that re- 
late quantum and classical correlation functions in terms of 
fluctuations of the density of states and correlations of par- 
ticle density are also examined and support the experimen- 
tal results. The results establish a correspondence between 
quantum spectral correlations and classical decay modes in 
an open systems. 

Pacs ref: 05.45.Mt, 05.45.Ac, 03.65.Sq, 84.40.-x 



1 Introduction 

The correspondence between quantum and classical me- 
chanics has turned out to be particularly rich for sys- 
tems that are classically chaotic fit]. This correspon- 
dence is usually considered in terms of the asymptotic 
behavior of quantum properties. Statistical analysis of 
the quantum spectra of closed chaotic systems display 
short-range correlations that are universal and can be 
described by random matrix theory (RMT) [||, ||| . The 
long range behavior of spectral fluctuations require semi- 
classical treatments that incorporate periodic orbits [Q. 

In this paper we establish a new connection between 
the statistics of quantum spectra and classical dynamical 
properties of an open chaotic system, utilizing physical 
and numerical experiments on the well-known n-disk bil- 
liard geometry. The n-disk system || is a paradigm of 
open hyperbolic systems, with an associated strange re- 
peller [[| whose strange properties arise from the Cantor 
set of the unstable periodic orbits (PO). For hyperbolic 
repellers, Ruelle and Pollicott || showed that the 
dynamical evolution of repellers toward equilibrium can 
be expressed in terms of complex poles, leading to the 
so-called Ruelle-Pollicott (RP) resonances. The n-disk 
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billiards have been extensively addressed theoretically, 
both classically and semiclassically || ||, |To[ Oh and 
only recently have been studied experimentally |l4 |l3ll . 

Our microwave experiments on the n-disk billiards 
measure the quantum transmission spectrum 6*21 (fc), 
which is essentially the two-point Green's function 
G(ri,r 2 ,fc), and from which the spectral autocorrela- 
tion C(n) is determined. The small wave vector k (long 
time) behavior of the spectral autocorrelation provides 
a measure of the quantum escape rate, and is in good 
agreement with the corresponding classical escape rate 
7c/ [|l2| as predicted by RMT analysis of this uni- 
versal behavior. For large k > 7 C ; (short time), non- 
universal oscillations of the autocorrelation C (n) are ob- 
served, that can be understood completely in terms of 
classical RP resonances. This behavior is similar to that 
of the variance S 2 (L)of closed chaotic systems, which 
also displays leading universal behavior for small L, fol- 
lowed by non-universal behavior for large L that can be 
described in terms of system specific periodic orbits Qj. 
The present work demonstrates that C{n) of open sys- 
tems displays a similar trend from short-range univer- 
sality to long-range non-universality, however the non- 
universal contributions to quantum correlations are de- 
scribed in terms of other classical phase space structures 
of the associated repeller, viz. the RP resonances [Q. 

A concise account of the main results was presented 
in Rcf. j^]. This paper further explores the details 
of the experimental aspects, as well as confirming the 
results via numerical experiments using simulated quan- 
tum spectra. Semiclassical arguments that demonstrate 
this correspondence are discussed in terms of fluctuations 
of the density of states and correlations of the particle 
density. The wider implications of the present results 
are also discussed. 



2 Experimental results 

The quantum dynamics of the 2-D n-disk system can 
be realized in a microwave experiment exploiting the 
mapping between the Helmholtz equation and the 
Schrodinger equation in their stationary forms and in the 
2-D limit. This is because under the conditions of the 
experiment, the Maxwell-Helmholtz equations reduce to 
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(V 2 + fc 2 )* = with 4" = E z the microwave electric 
field. The experimental geometry consists of thin cop- 
per disks sandwiched between two large copper plates of 
size 55 x 55 cm 2 in area. In order to simulate an infi- 
nite system, microwave absorber material ECCOSORB 
AN- 77 was placed between the plates at the edges. Mi- 
crowaves were coupled in and out through antennas in- 
serted in the vicinity of the scatterers. All measurements 
were carried out using an HP8510B vector network an- 
alyzer which measured the complex transmission ampli- 
tude (52i ) and reflection amplitude (5n) 5-parameters 
of the coax + scatterer system. An example of the mea- 
sured transmission T(k) — |5 2 i(fc)| 2 of a 4-disk system 
is shown in Fig. 1. For more details of the microwave 
experiments, see Ref. 

An open system can be represented by an effective 
Hamiltonian consisting of two parts, H = H c + iW, 
where H c is the Hamiltonian for the closed system 
and W — represents the decay to open chan- 

nels. Since the total Hamiltonian is not Hermitian, 
the eigenfunctions ip n of H with eigenvalues E n do not 
form an orthogonal set. Instead, let <j) n be the eigen- 
function of the adjoint operator H* with eigenvalues 
.E*, and the ip's and </>'s form a bi-orthonormal set 
with J ipm(r)(f> n (r)dr = S mn . The Green's function for 

an open system is flfj G(fi,r 2 ,k) = Y, n e-'e • 
For the n-disk system, with the quantum resonances in 
wave vector space k n = s n + is' n , we have 



\G(n,f 2 ,k)f 



\ipn(n)Mr2)\ 2 



2m \ -> 

~tf^[(k + s n y + s>Mk-s n y + s'*] 

h 2 ^ ' 

n.m 



The second sum with the prime includes the off-diagonal 
terms with n ^ m. Their contribution can be neglected 
if all the resonances are well separated. 

The trans- 
mission amplitude 52i(fc) measured between two point 
antennas is essentially the Green's function G(ri,r2,fc) 
of the system, so that 52i(fc) = A(k)G(fi, r 2 , k) with 
A{k) a slowly- varying modulating function Jl3fl . Setting 
c„ = (2m/h 2 )\A(k)ip n (f 1 )^ n (f 2 )\ 2 /[(k + s n y r +s% and 
ignoring the off-diagonal terms, we get the transmission 
coefficient T(k) expressed as a sum of Lorentzian peaks 



T(k) 



(k - s n ) 2 + s% ' 



(1) 



The coupling c n in Eq. (Q) depend on the location of 
the probes. The k dependence of c n is locally weak 
near a resonance since l/[(k + s n ) 2 + s' 2 } is flat around 
I: ~ s„ » s' n for sharp resonances. One can treat c n 
practically as constants in the neighborhood of a reso- 
nance. This shows that the measured transmission co- 



efficient T(k) shown in Fig. 1 directly yields the wave- 
vector s n and the half-width s' n of the sharper scattering 
resonances. We next compare the experimental reso- 
nances with semiclassical calculations. 

2.1 Calculation of Quantum Resonances 

Although any chaotic system can be quantized numer- 
ically by directly diagonalizing the Hamiltonian trun- 
cated in certain expansion space, for a large number of 
systems, one uses the techniques based upon semiclas- 
sical periodic orbit theory such as the cycle expan- 
sion frij| , Fredholm determinant |]l8f or harmonic inver- 
sion Except for the lowest eigenvalues, the semi- 
classical quantization gives very accurate results. Using 
Gutzwiller's trace formula JjJ, the semiclassical Ruelle 
zeta-function can be derived as an Euler product over 
all the prime PO which are the PO without repetition 

@0 



Cj,sc — J |_ ( ^ tp.se) 1 



with t PjSC the semiclassical weight 



t p , sc = (-!)'» exppfcLpVIApl 1 / 2 ^. 



(2) 



(3) 



Here l p is the number of collisions of the PO with the 
disks, L p is the length of the PO and A p the eigenvalue 
of the instability matrix J p . j comes in from the decom- 
position of the trace of the Green's function into Ru- 
elle zeta- functions with j = 0, 1, • • •. The quantum reso- 
nances are the poles of the Ruelle zeta-functions which 
is directly related to the trace of scattering matrix or the 
trace of the Green's function. If a symbolic dynamics ex- 
ists for the system, only a few prime POs will be needed 
in the zeta-function to give quite accurate eigenvalues 
because the curvature in the cycle expansion will decay 
exponentially with increasing PO length [ p~7| . 

As an example, consider the non-chaotic 2-disk sys- 
tem with disk separation R and disk radius a. There 
is only one unstable prime PO between the disks. The 
symmetry group of PO is C 2 , with two one-dimensional 
irreducible representations, symmetric A\ and anti- 
symmetric A 2 . The semiclassical resonances in wave 
vector space are k n = [nir + i(l/2) In A]/(R — 2a) with 
A = a — 1 + -Jcr(a — 2) the eigenvalue of the instability 
matrix in the fundamental domain and a = R/a. Here 
n is odd for A\ representation, n even and n ^ for A 2 
representation p0[ . 

For the chaotic n-disk system with n > 3, there is no 
analytical expression of the semiclassical quantum reso- 
nances. To do the numerical calculations of the quantum 
resonances, use can be made of the symmetry of the sys- 
tem and the cycle expansion in the fundamental domain 
[[To[ 11 . We have calculated numerically the quantum 
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resonances of the 3-disk and 4-disk systems with differ- 
ent disk separation ratio a by searching the poles of the 
Ruelle zeta function (||). An example of the results of 
the calculation for the s n , s' n for a 3-disk system is shown 
in Fig. 2. The experiments only access a small range of 
the quantum resonances with < Re(/) < 20GHz and 
—0.3GHz < Im(/) < 0. The experimental quantum 
resonances are well matched by the sharp resonances 
from the semiclassical calculations. The very broad res- 
onances are not observed because the coupling to them 
is weak. For details concerning the comparison of ex- 
perimental and calculated quantum resonances, see Ref. 

2.2 Experimental Spectral Autocorrela- 
tion C(k) 

The spectral autocorrelation function is determined from 
the experimental 
spectra as C t (k) = (T(k - (k/2))T(A; + (/c/2))) fc with 
average carried out over a band of wave vector centered 
at certain value ko and of width Ak jL3]. The spectral 
autocorrelation can be fitted well by a superposition of 
Lorentzians : 



C t {k)= £ 



- i7l '2 + ( K±7 ;0 2 



(4) 



with Y n ± i 7 ^' the classical RP resonances in wave vector 
space. A semiclassical derivation will be given in Sec. 3. 

The experimental RP resonances are obtained by fit- 
ting the experimental autocorrelation as shown in Fig. 3 
with Eq. (^). Similar comparisons between experiment 
and theory for the n-disk system with n = 2,3,4 are 
presented in Ref. |l5). The coupling b„ in the above 
equation which should depend on the location of the 
probes determine the decay probability of the classical 
RP resonances. Since we do not have knowledge of them 
experimentally, they are chosen to optimize the fitting. 

The experimental RP resonances (7^,7") obtained 
from the Lorentzian decomposition of the experimen- 
tal C(k) as shown in Fig. || are displayed in Table 1. 
They can be compared with theoretically calculated RP 
resonances, which are the poles of the classical Ruelle 
zeta-function, and whose calculation will be explained in 
the next subsection. Although the position of the peaks 
of the oscillations in the experimental autocorrelation in 
Fig. |3] is quite accurately given by the imaginary part 
7^' of the RP resonances, the experimental half-width of 
the oscillation is almost always smaller than the calcu- 
lated real part 7^ of the RP resonances. This is possi- 
bly because the absorber used is not ideal, and leads to 
quantum resonances that are slightly sharper than cal- 
culated. This results in systematically sharper widths of 
the experimental RP resonances. 



2.3 Calculation of Classical resonances 

The theoretical RP resonances are essential ingredients 
of a Liouvillian description of the classical dynamics, and 
can be calculated as the poles of the classical Ruelle zeta- 
function. The classical dynamics of Hamiltonian flow in 
symplectic phase space can be described by the action 
of a linear evolution operator which is also the Perron- 
Frobenius operator £'(?/, x) = 5(y — f t (x)). Here /'(a?) is 
the trajectory of the initial point x in phase space. The 
trace of the Perron- Frobenius operator is tr£* = J 6(x — 

f(x))dx = j: p T p S(t - rT p ) |det(l - Jp|~\ The 
first summation is over all the prime classical PO p with 
period T p , the second one is over the repetition r of the 
prime PO. 3 p again is the instability matrix. The evolu- 
tion operator has a Lie group structure £* = e~ At with 
A the generator of the Hamiltonian flow. This gener- 
ator A has the classical resonances as its eigenvalues, 
7 „ = 7,'j ± 17", which we call the RP resonances 0, ||]. 
We have tr£* = J2^Lo ff™ e_7 "*i with g n the multiplicity 
of the resonances. These resonances determine the time 
evolution of any classical quantity. The trace formula for 
classical flows is obtained from the Laplace transform of 
the above expression |pl| 

/>oc 

tr£(s) = / cfte st tr£* = tr(A - si)' 1 
Jo 



~i |det(l-J;) 



(•5) 



In order to calculate the RP resonances, the Ruelle 
zeta-function is introduced as 



Cf3,cl — — tp,cl) 1 ■ 



(6) 



Here the product is over all the prime PO with t PjC i the 
classical weights of the periodic orbits 



t p , cl = exptsTpVIAplA^" 1 



(7) 



Here the integer (3 > comes in from the expansion 
of the determinant |det(l — J£)| . One can see that 
the above classical Ruelle zeta-function is very similar 
to the semiclassical one (Eq.(g)). The classical Ruelle 
zeta-function is exact and can be directly derived from 
the above trace formula (|5|) of the Perron-Frobenius op- 
erator. The topological pressure P((3) can be defined 
as the simple pole of the classical Ruelle zeta-function 
[ p2[ with P extrapolated to the entire real space. All 
the characteristic quantities of the classical dynamics, 
such as the Kolmogorov-Sinai entropy, escape rate, frac- 
tal dimensions, can be derived from P((3) j|. For the 
hard-disk system, the classical velocity v of particles is 
constant. So T p = L p /v. For simplicity, one can set 
v = 1. The RP resonances are then calculated in wave- 
vector space. 
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For example, the classical RP resonances of the 2-disk 
system are 7„ = [In A+inir]/ (R—2a) with A given before. 
Here n is even for A\ representation and n odd for A 2 
representation. Similar to the situation of the quantum 
resonances s n + is' n , there is no analytical expression of 
the RP resonances 7^ ±27^' for the chaotic n-disk system. 
They can be calculated numerically as the poles of the 
classical Ruelle zeta- function Eq. (||) ||. 14 prime PO 
up to period 3 were used in our calculation of the RP 
resonances of the 4-disk system, the results of which are 
shown in Table 1. 



2.4 Numerical simulation 

To gain further insight into the established experimen- 
tal correspondence, we have utilized the results of the 
calculations of the quantum and classical resonances to 
carry out a numerical simulation of the physical exper- 
iment. We calculate numerically a close approximation 
to the experimental quantum spectrum using the quan- 
tum resonances. We them determine the spectral auto- 
correlation of the simulated spectrum, and compare its 
decomposition with the calculated classical resonances. 

Writing Eq. (|l|) explicitly, we have following form of 
the transmission T(k) ~ 

(2m/h 2 )\A(k)\ 2 £„ m .$4%$%+^ ■ Here we 
assume the presence of the probes poses a small per- 
turbation and will not change the quantum spectrum. 
For a closed chaotic system with time-reversal symme- 
try, the wave density is expected to follow the Porter- 
Thomas (PT) distribution in RMT (23). For a very 
open system and when the probes are far away from 
the scatterers, the wave density will not follow the PT 
distribution. But in the vicinity of the scatterers, one 
may assume the density to follow the PT distribution. 
Setting pi n = \ip n (ri)\ 2 and p 2n = |0™(r2)| 2 , one has, 
P(pn) = (2np n p) exp(-p„/2p) with p being the 
ensemble averaged density. The densities at the loca- 
tion of different antennas were found not to be cor- 
related in closed cavities 



(2m/h 2 )p ln p 2n = 1. We thus get 



T(k) 



24 . We assume the same 



is true for an open chaotic system in the vicinity of 
the scatterer disks. For simplicity, we consider the en- 
semble average of the transmission coefficient T(k) ~ 
(2m/h 2 )\A(k)\ 2 p 1 p 2 £„ 1/ [(k + s n ) 2 + s' 2 ] [(k - s n ) 2 + 
From this expression, one can see that the contribution 
from broad resonances is suppressed. The main con- 
tribution is from sharp resonances. This has already 
been observed in our experiments p^, For these 

sharp resonances, if they are far away from the ori- 
gin, (k + s„) 2 + s' 2 ~ 4fc 2 . Thus, for large k > s ' n , 
the above transmission can be approximated as T(k) ~ 
(2m/h 2 )\A(k)\ 2 (p 1 p 2 /4k 2 ) E„ [(* ~ s n ) 2 + s' 2 ]-\ The 
function A(k) is found to be proportional to k |25| , 
p6f. For convenience, we set A(k) — 2k and also 



[{k - s n ) 2 + s>- 



(8) 



The k dependence of A(k) can be understand as follows. 
In the experimental setup, one of the antenna can be re- 
garded as a diploe radiating electromagnetic waves. The 
radiated electrical field of a dipole is proportional to k if 
an alternating current with constant amplitude I- m was 
maintained on it |2^| . The voltage picked up on another 
antenna is just the electrical field at the antenna location 
times the length d of the antenna inside the cavity. The 
transmission amplitude S 2 \ measured by the analyzier is 
the ratio of the output voltage on one antenna to the in- 
put current of another antenna, S 2 \ = dE ZtOUt / Ii n , thus 
A{k) = dE zM /I in oc k. 

The expression (|§|) is thus used to get the averaged 
transmission and then to calculate the autocorrelation 
Ct(k) numerically. We have done that for the n-disk 
system with n = 2,3,4 using a = 5cm. Here we 
just present the results for the 3-disk system as shown 
in Fig. 4 and 5. Note that / = {c/2n)k, 1GHz is 
equivalent to 0.2096cm . The symmetry group for 
the chaotic 3-disk system is C^ v JlQ] , which has two 
one-dimensional irreducible representations A\, A 2 , and 
one two-dimensional irreducible representation E. In 
the fundamental domain, only the resonances of the A 2 
representation will contribute to the transmission |l2| . 
The semiclassical resonances of A 2 representation in the 
range < Keka < 100 and —0.5 < Imka < for the 
disk separation R/a — 4\/3 are obtained from Q c (ik) 
with j — 0. About 200 resonances are obtained as shown 
in Fig. 2. 

The corresponding classical resonances shown in Fig. 
5 are from the Ai representation of the classical Ruelle 
zeta- function (g) with (3=1. 8 prime PO up to pe- 
riod 4 were used in our calculations of the quantum and 
RP resonances. The numerical transmission and the au- 
tocorrelation are shown in Fig. 4 and 5, respectively. 
One can see that the oscillations in the autocorrelation 
are fully determined by the classical RP resonances as 
predicted by Eq.([|). 

,<i. Quantum and classical correla- 
tions 

In the above section, we have demonstrated experimen- 
tally and numerically that the quantum auto-correlation 
and hence the statistics of quantum resonances are deter- 
mined by the RP resonances of the underlying classical 
dynamics. Here we explore the theoretical justification 
for this correspondence. 

For chaotic open scattering systems, the distribution 
and correlation of the quantum resonances have been 
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discussed in the framework of RMT 28 . However, the 
RMT results are the ensemble averages of different sys- 
tems with certain symmetry in common. Consequently, 
RMT does not make contact with system specific fea- 
tures of the underlying classical dynamics. Thus it is un- 
able to make any concrete connection between the quan- 
tum and classical dynamics, for which it is necessary to 
go beyond RMT. 

The appearance of nonuniversality requires the use of 
structures such as periodic orbits within semiclassical 
calculations. To proceed, we first consider the correla- 
tion of the density of states. In the semiclassical periodic 
orbit theory, the fluctuation part of the density of states 
is 



V{k)- 



■ Im tr£r(ri, r%, k) 



irkL v 



det(l - J? 



1 1/2 ' 



(9) 



The trace of the Green's function is the integration over 
the configuration space with r\ = r*2. The autocorrela- 
tion is 



1 

2K 



K 



K 



dkV(k + ^)V(k- -), K-»k. 



Substituting the semiclassical expression d) ofD(fc) into 
the above equation, one gets 



\ |det(l- j;) 
in 



Retr 



A — si 



(10) 



with s = ivk. Since the spectrum of the generator A is 
u(7^ ± ij'n), one can see that in the semiclassical theory, 
the autocorrelation of the density of states is determined 
solely by the classical RP resonances. The transmission 
spectrum T(k) is a projection of the density of states, 
one expects that the autocorrelation of the transmission 
is also determined by the RP resonances. The approxi- 
mation in Eq. ( |To| ) is valid only for hyperbolic systems 
with finite symbolic dynamics. For these systems as ex- 
plained in Sec. 2.1, the contribution to the autocorre- 
lation Cd(k) is mainly from a few fundamental prime 
POs. For those PO, one has T p 2 ~ TT p with T the av- 
erage period. For weakly open or closed system, the 
approximation will break down. 

The experimental quantity whose correlations we ex- 
amined in Sec. 2 is the Green's function, while theo- 
retical arguments usually examine the density of states 
or the time delay |]30f which are however not directly 
measurable experimentally. To develop the correlation 



theory for the transmission we measured, consider the 
stationary Green's function G(f 7 r ,k) which is the so- 
lution of the following equation (V 2 + k 2 )G(r, r , k) = 
S(f — ro), with certain boundary condition. The 
quantum mechanical propagator is the Fourier trans- 
form of the Green's function. One can construct 
a time evolution of a wave packet K^fo ,t) with 

k(f,f Q ,t) = (2-kM)- 1 f Ae G(f,r , yj2me/n 2 )e- t£t ^de, 

with e = H 2 k 2 /2m. Here the integration is per- 
formed around eo = h 2 k 2 /2m, with the range Ae = 
HvAk and v = Kko/m the group velocity of the 
wave packet. The integration in the e space can be 
changed into that in the k space. We get K(f,fo,t) = 
e- t£ot / h (v/2iri) J Ak G(r, f ,k + k)&- mkt dk. The propa- 
gator K(j, ro, t) is just the wave function at point rdue 
to a (5-function excitation of the system at point r*o and 
time t = 0. The particle density p(t) is thus 



K(r,f ,t) = -—t / G(r,f ,k 
G*{f,rQ,k + k')e mk,t dk', 



k)e 



{) ki 



dk 



Pit)- 



Ak 

with 

the average given by (p(t)) t = lim-r^oo T 
{v 2 /2ttL)J dk\G(k)\ 2 . A probability measure is as- 
sumed for the above average to converge 0. Here L 
is the size of the whole system. For an open system, 
L — ► oo. Here we used the definition of the S function 
S(t) = (27T)- 1 J due-*"* . 

The autocorrelation of the particle density is 



JoP(t)dt 



C p (r) = (p(t)p(t + r)) t 



(P)l 



Using the diagonal approximation, we get 



C P {r) = 



4tt 2 L 2 



c?kCt(k)6 



(11) 



Note that we have T(k) oc |G(r, r , k)\ 2 . 

If one assumes that the correlation of the particle den- 
sity is classical, one has [M] 



C P {r) 



= ^26;e^> T cos 7 ;W, 



(12) 



where the coefficients b' n are the coupling, and j' n ± 27" 
the RP resonances evaluated in wave- vector space. The 
Fourier transformation of ([ll]) and (O) gives 



b ni n 



(13) 



with b n = 2irL 2 b' n /v' l Ak. The autocorrelation is usually 
normalized as Gt(0) = 1 so that the system size L will 



5 



be cancelled. Thus the above expression is valid for both 
closed and open systems. 

The above arguments closely parallel Agam's |}2| re- 
cent application of the diagrammatic techniques for dis- 
ordered systems to open chaotic systems. We remark 
that the above arguments differ somewhat from those of 
Agam, whose results are valid when the two probes are 
far away from each other which implies {T 2 ) = 2 (T) 2 . 
This is the case in the diffusive limit with a Poisson dis- 
tribution of the transmission coefficient T. 

The above correspondence between the autocorrela- 
tion derived from semiclassical theory and the classical 
resonances is obviously not exact, but requires correc- 
tions due to quantum interference. These corrections are 
small for open systems, but are expected to dominate for 
closed systems. 



4 Discussion and concluding re- 
marks 

The results presented here establish an alternative path 
from quantum to classical achieved by considering the 
correspondence of quantum correlations and classical de- 
cay modes. As Fig. || shows, the small K behavior of 
C(k) is universal, while the large k behavior shows non- 
universal oscillations that are completely described by 
the classical RP resonances. This demonstrates that in 
an open system, quantum correlation functions are inti- 
mately related to decay laws of classical survival proba- 
bilities or evolution of particle ensembles. 

Most theoretical calculations in quantum chaos us- 
ing ensemble averaging are restricted only to the lead- 
ing order consistent with RMT, which probes only the 
average properties of the system, such as the aver- 
age escape rate 7. The leading RMT contribution to 
the correlation C(k) is Lorentzian ||| |34|, 35, [l^] 
C(k) = C(0)/(1 + (k/j c i) 2 ), when the classical decay is 



exponential, which occurs for hyperbolic systems such as 
the n-disk billiards. Our earlier experiments established 
this result in the open microwave billiards jl^] . 

As we have shown, our experiments go well beyond the 
average information represented by RMT - the spectral 
data contain further details of the fine structure embod- 
ied in the decay modes or classical resonances of the sys- 
tem. In terms of the fractal repeller, the wave mechan- 
ical experiments thus yield not only the leading average 
property 7 C ; related to the fractal dimension, but also 
further higher order properties of the repeller in terms 
of the classical resonances. While we have demonstrated 
this for the n-disk billiard, the above result should be 
valid generally for hyperbolic open chaotic systems. 

The relevant classical limit here is the evolution of 
classical particle ensembles. The notion of ensemble av- 
eraging that is common in quantum mechanics is also 



useful in classical treatments of chaotic systems. This 
is because while the prediction of individual particles is 
sensitive to initial conditions, the dynamics of appro- 
priate averages of a perturbed system relaxing towards 
equilibrium is well defined. The RP resonances deter- 
mine the decay modes of the dynamical systems evolv- 
ing from a nonequilibrium state to equilibrium or steady 
state. For an open system, Gaspard and Nicolis |3^] 
have derived an escape rate formula to relate the co- 
efficient of diffusion and the escape rate of large open 
systems. While for an open system, the escape rate 70 
is nonvanishing, for closed system, 70 = 0, and the cor- 
responding eigenstate is just the equilibrium state. The 
RP resonances are essential ingredients of a description 
of irreversibility and thermodynamics based upon mi- 
croscopic chaos j!o[ ||, [IT]. However these resonances 
have remained mathematical objects. The present ex- 
periments have provided physical reality to them. It is 
remarkable that we are able to observe classical reso- 
nances in a quantum experiment. Because microwaves 
do not interact (collide) with each other, ensemble av- 
eraging is easier to achieve in microwaves than in fluids. 
Thus theories that employ ensemble averaging are ide- 
ally applicable to the microwave experiments. This is 
true for both open systems as shown here, and for dis- 
ordered billiards |Q. 

The present work shows that the information concern- 
ing the classical resonances is coded into the quantum 
spectrum. This is somewhat anticipated in periodic or- 
bit theory, where one can observe the strong similarity 
between the classical (Eq.(^|) and ([?])) and semiclassi- 
cal (Eq.(||) and (||)) Ruelle zeta- functions. A one-to-one 
correspondence exists between k n and 7„ for the 2-disk 
system as discussed in Sec. I. Another example where an 
exact relation between the classical and quantum reso- 
nances exists is the system of a free particle sliding on a 
compact surface of constant negative curvature, the RP 
resonances were found by Biswas and Sinha [j37f to be 
70 = 0, 7„ = \ ± ip n for n > 1 with real p n the solution 
of a certain equation while the quantal excitation eigen- 
values are E n — p^ + j. One would also expect a relation 
between the distribution of the classical and semiclassi- 
cal resonances. The RP resonances have recently been 
observed in the distribution of the zeroes of the Riemann 
zeta function |f38f . 

Our work demonstrates that suitable quantum corre- 
lations diffuse just like classical observables in an open 
system. An interesting question is what will happen 
when the number n of disks are very large, a situation 
that is frequently referred to as the Lorentz gas. Due to 
the exponential separation of the nearby orbits, a par- 
ticle is forced to move through the bouncing paths, or 
undergo diffusive motion, for long time. One may argue 
that the particle motion occurs in a system of size L, 
in which it has a diffusive path and hence the average 
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escape time t c ~ L 2 /D, with D the diffusion coefficient. 
Thus the escape rate 7 ~ D /L 2 . A more general expres- 
sion for diffusion coefficient D is derived through the es- 
cape rate formula in Ref . [^9) . For the Lorentz gas with 
large L, by taking the general form of the escape rate 
7 = A — hks , where A is the Lyapunov exponent and h^s 
is the Kolmogorov-Sinai entropy of the fractal repeller 
<F . l > one has D = [L/2.40482] 2 [\{T L ) - h KS {T L )]. 

There is a fundamental connection between the 
present observation of classical resonances and the dif- 
fusion operator approach of the supersymmetry theories 
[ [f3t of disordered systems. The statistics of the spec- 
tra and eigenfunction of disordered systems is directly 
related to the classical spectral determinant D(e) = 
e~ 2 11^(1 + e2 ^ 2 /7^) _1 : where e is the energy ( in unit 
of A, the average level spacing), and j u are the non-zero 
eigenvalues of the diffusion operator Q, |i| . These are 
just the RP resonances that we have observed in open 
systems. Using the supersymmetric techniques, Agam, 
Altshuler, and Andreev jl6) calculated the level correla- 
tion and showed that it can be fully determined by the 
classical spectrum of the Perron-Frobenius operator. Bo- 
gomolny and Keating p7| arrived at a similar conclusion 
using the periodic orbit approach. The applicability of 
these theories for closed systems is currently a matter of 
debate |Q, but the applicability of the supersymmetry 
theory to disordered systems is well established, and is 
strongly supported by microwave experiments on disor- 
dered billiards Q ff|] . 

As pointed out by Berry the two limits T — > 00 
and h — > do not commute with each other, thus 
the long-time quantum evolution is fundamentally dif- 
ferent from long time classical evolution. The proba- 
bility P q that the quantum chaotic system decays at 
time t after its formation follows a power law |24| 
P q (t) ~ (1 + 2Tt)~ 2 ~~ , for t larger than the equilibrium 
time of the system. Here T is the average decay width, 
M the number of open channels. As the system is opened 
up with more channels, the decay will be almost expo- 
nential since the above equation can be approximated by 
P q {t) ~ exp [-(4 + M)Tt] , for large M and the weight 
of the algebraic tail becomes negligible. More precisely, 
as pointed out in Ref. [pO] , there exists a new quantum 
relaxation time scale t q <tu, with tn = h/A the Heisen- 
berg time. Beyond that time, the quantum correlation 
decays algebraically. For shorter times, the quantum 
correlation decays exponentially in the same way as the 
classical correlation does. For an open system, since the 
spectrum is continuous, the Heisenberg time is actually 
infinite, and so is the quantum relaxation time t q . This 
is one reason that the classical resonances are clearly 
visible in the present experiments, and the Lorentzian 
decomposition of the spectral autocorrelation works so 
well. 

The transmission spectrum measured in the mi- 



crowave experiment directly corresponds to the conduc- 
tivity measured in electronic quantum dots. In quantum 
dot experiments though, the magnetoconductivity a(B) 
is typically measured, so that the resulting autocorrela- 
tion C(AB) is related to the area spectrum, rather than 
the wave vector dependence C(/c) measured here which is 
related to the length spectrum [Q . Nevertheless similar 
arguments should be applicable there also, as well as to 
a variety of other physical systems which are essentially 
open and chaotic, such as molecular photo-dissociation 
and chemical reactions || . 
While the coupling to the quantum resonances is well 
understood, the statistics and also the physical mean- 
ing of the experimental coupling b n to the classical reso- 
nances in Eq. ([!]) deserves further clarification. It would 
also be very interesting to examine the role of classical 



resonances in mixed phase space systems 52 , such as 
the open billiards frequently considered in quantum dots 
where the classical decay is non-exponential |53 . More 



detailed work should be devoted to these and other in- 
teresting questions raised by the present results. 
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Table 1: Experimental and theoretical RP resonances of the 
4-disk system in the fundamental domain with R/a — 4%/2 
and a = 5 cm. The RP resonances are in the unit of cm - . 
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Figure 2: Quantum resonances in the A2 representation of 
the 3-disk system with R/a — 4^3 with a = 5 cm in the 
fundamental domain. The quantum resonances we observed 
experimentally are in the range < / < 20GHz. See for 
details. 
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Figure 1: Experimental transmission T(k) of the 4-disk 
system in the fundamental domain with R/a = 4\/2 and 
a — 5 cm. Here / = (c/2n)k. Inset: Geometry of the 3- 
and 4-disk billard. The solid lines represent the fundamental 
domain in which the experiments were carried out. 
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Figure 3: Autocorrelation function C(k) vs k (cm -1 ) of the 
4-disk system in the fundamental domain with R/a = 4%/2 
and a — 5 cm. The gray line is the experimental C(k) cal- 
culated from the experimental trace shown in Fig. 1. The 
data show the small k universal behavior followed by the 
non-universal oscillations for large k. The entire autocorre- 
lation C(k) can be described in terms of the RP resonances 
using Eq.(4) (thick solid line). The thin lines represent the 
Lorentzian decomposition into individual RP resonances. 
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Figure 4: Numerical spectral transmission T(k) in the fun- 
damental domain of the 3-disk system for R/a = 4\/3 and 
a — 5 cm. 




Figure 5: Numerical spectral autocorrelation C(n) (top) 
and RP resonances 7' ± ry" (bottom) for the 3-disk system 
with R/a — 4\/3 and a — 5cm in the fundamental domain. 
C(k) was calculated from the numerical spectrum shown in 
Fig. 4 for interval Ak — 6cm" 1 with the central values ko in 
cm -1 indicated in the figure. 
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